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We consider a two-leg Bose-Hubbard ladder in the presence of a magnetic flux. We make use of 
Gross-Pitaevskii, Bogoliubov, bosonization, and renormalization group approaches to reveal a struc¬ 
ture of ground-state phase diagrams in a weak-coupling regime relevant to cold atom experiments. 

It is found that except for a certain flux (j) = tt, the system shows different properties as chang¬ 
ing hoppings, which also leads to a quantum phase transition similar to the ferromagnetic XXZ 
model. This implies that population-imbalance instability occurs for certain parameter regimes. On 
the other hand, for cf = n, it is shown that an umklapp process caused by commensurability of a 
magnetic flux stabilizes a superfluid with chirality and the system does not experience such a phase 
transition. 

PACS numbers: 67.85.-d,05.30Jp 


I. INTRODUCTION 

Quantum systems subject to high magnetic heids are 
known to acquire nontriviai characteristics such as the 
Hofstadter butterfly [T] and quantum Haii effect [2] in two 
dimensionai systems. Recentiy, the so-caiied synthetic 
gauge fieids [3H5] avaiiabie in uitracoid atomic systems 
pave the way to reaiizations of such systems. An advan¬ 
tage of coid atoms is that one can controi the geometry, 
dimension or quantum statistics of systems and param¬ 
eters of microscopic Hamiitonians at unprecedented iev- 
eis, which is used to expiore non-triviai quantum states. 
Aiong these iines of researches, the Hofstadter Hamiito- 
nian 011 ] and Haidane topoiogicai modei 0 have been 
reaiized in coid atoms. 

In addition to such a non-triviai feature of non¬ 
interacting quantum matters, as is weii known, an in¬ 
teraction is a key ingredient for the diversity of nature. 
Indeed various phenomena such as superconductivity, su¬ 
perfluidity and Mott transition are understood as a con¬ 
sequence of the interactions. One can thus expect the 
interaction effects in such topological matters to cause 
further non-triviai nature, and effective approaches in¬ 
corporating interactions are required in theory. 

In one dimension one can successfully apply field theo¬ 
retical approaches to incorporate interaction effects in a 
non-perturbative manner. Thus, reduction of dimensions 
would be a way to understand physics involving magnetic 
fields and correlations. The minimal model to show non¬ 
trivial effects in the presence of magnetic fields is a two- 
leg ladder, and the bosonic version was first discussed 
in Ref. 0 in the context of Josephson junction arrays. 
In this study, it has been predicted that two different 
phases show up: Meissner and vortex phases. While in 
the former phase a chiral current analogous to a Meissner 
edge current is induced on the legs by a magnetic flux, 
in the latter it is significantly reduced due to penetration 
of vortices, which is analogous to field-induced vortices 
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FIG. 1. A schematic figure of a two-leg Bose-Hubbard model 
with a flux (j). In the gauge chosen here, the flux effect appears 
only in rung hoppings. 


in type-II superconductors. Later on, the theoretical in¬ 
terest has been devoted to the strong correlated regime, 
and it has been demonstrated that commensurability of 
a particle filling poses a Mott insulator with chirality and 
interesting critical properties HMD- 

Remarkably, the two-leg bosonic ladder subject to a 
magnetic flux has been successfully realized in an ex¬ 
periment m, where a weakly-interacting regime is con¬ 
cerned. In this experiment, it has been confirmed that 
behaviors of chiral currents are consistent with what has 
been predicted in Ref. 0. Thus it seems that a basic 
consensus in a weak-coupling regime is obtained. 

More recently, however, it has been argued in Ref. [T^ 
that in a weak-coupling regime there should exist an ad¬ 
ditional phase where a spontaneous population imbalance 
between the legs occurs. This additional phase named a 
biased ladder phase has been shown with the theory of 
weakly-interacting Bose gases normally used in higher di¬ 
mensions m At the same time, in one dimension quan¬ 
tum fluctuations should be non-negligible in most cases. 
Thus it is worth considering whether the biased ladder 
phase is still robust against quantum fluctuations. 

In this paper, we examine the two-leg Bose-Hubbard 
model in the presence of a magnetic flux in a weakly- 
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interacting regime by means of a couple of effective- 
theory approaches. We show that a spontaneous pop¬ 
ulation imbalance indeed occurs, and is stable against 
quantum fluctuation effects. In particular, we point out 
that the effective theory has the similarity to a ferro¬ 
magnetic XXZ quantum spin model. This implies that 
the Heisenberg point exists in the phase diagram, where 
SU{2) symmetry shows up in the low-energy effective 
theory although the original Hamiltonian does not pos¬ 
sess that symmetry. This situation is somewhat similar 
to a two-leg extended Bose-Hubbard system analyzed in 
Ref [IH], where the low-energy effective theory possesses 
an emergent symmetry. 

We also state that umklapp processes coming from 
commensurability of a flux should be seriously consid¬ 
ered at the mean-field level, which has been overlooked 
in the previous studies. The umklapp processes exist¬ 
ing in (/) = TT destabilize the biased ladder phase, and as 
a consequence, only the commensurate vortex phase is 
allowed. 

The structure of the paper is as follows. In Sec. [IT] we 
review structures of single-particle bands as a function 
of 4>, and low-energy effective Hamiltonian reflecting the 
band structure. The single-particle band bottom shows 
different topologies: A single minimum for a small flux, 


and double minima for a large flux. In Secs. HI and IV 


we discuss physics separately for a single-minimum and 
for a double-minimum band structure, in which Meissner, 
vortex, and biased ladder phases are allowed depending 
on (j) and Kj J. In Sec. [V| summary and perspective on a 
phase transition between the Meissner and biased ladder 
phases, and on a stronger interaction effect are provided. 
Technical details on renormalization group equations are 
addressed in the Appendix. 


II. FORMULATION OF THE PROBLEM 

Following the setup in Ref. m, we consider the fol¬ 
lowing two-leg Bose-Hubbard ladder Hamiltonian: 

1=1 p=l,2 

+ H.c.) + - lOl) 

i i,p 

where J and K are the hopping amplitudes along the 
leg and rung directions, respectively. The applied flux 
is introduced via the Peierls substitution, and the corre¬ 
sponding gauge fields along the chain and rung directions 
are denoted by ^ and A^, respectively. The flux (j) is 
then given as (j) = ^ The tech¬ 

nology of laser-assisted tunneling [ilZllIMI] generates 
spatially-dependent phase in rung hoppings, which leads 
to the (j) flux per plaquette as described in Fig.fTl Namely, 
in this paper we choose the following gauge as A"^ = 0 


and Af- = (j)l. Taking into account the fact that an in¬ 
teratomic interaction is given by an s-wave scattering 
length, and the stability in bosonic systems, we restrict 
ourselves to a local repulsive interaction, U > 0. 

Apparently, the Hamiltonian Q is invariant under the 
simultaneous transformations, &i,i( 2 ) —>■ h, 2 {i) and 4> —>■ 
—(j). Thus, we can safely take the domain of definition in 
(j) as 0 < (j) < n. 

By using the general relation between current and 
Hamiltonian, j = — where A is the gauge field, we 
can define current operators along the legs and rungs as 

~ '^’^i^l+i,p^i,p ~ ^l,p^i+i,p)^ (2) 

jt = iK . (3) 

For the sake of convenience, we also introduce chiral cur¬ 
rent along the legs: 

= ( 4 ) 

We will see that and play important roles in char¬ 
acterizing each phase. 

Since we are interested in the regime J,K ^ U, 
we start with diagonalizing the single-particle Hamil¬ 
tonian. To this end, we perform gauge and Fourier 
transformations as bi^i = ^ h ,2 = 

^ Then, by considering a unitary 

transformation for hk,i and bk ^2 



the single-particle Hamiltonian can be diagonalized [HI 

IHj as 


Ho = '^{E+{k)alak +E_{k)l3l/3k). ( 8 ) 

k 

Here, the single-particle spectrum is given by 
E±{k) = 2J ~^ 4= Y -I-sin^ sin^ fc 


( 9 ) 
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which depicts the two-band structure as well as the 27r 
periodicity, reflecting the two-leg ladder geometry. 

Unless a strong interaction is concerned, single-particle 
low-energy states, bottoms in the lowest band, play im¬ 
portant roles in the low-energy many-body states. With 
this understanding, we neglect effects of the higher band 
Qffc by keeping in mind the condition J, K ^ U. 



FIG. 2. Changes of topology in the lower band at a certain 
(j> ^ TT. The direction of the arrows means the reduction of 
K/J. In a strong enough K/ J, the band has a single minimum 
while in the opposite limit, the band of double well structure 
forms. In between, there exists a critical point in which the 
band bottom becomes quartic in k. The double-well structure 
is always maintained at = tt regardless of values of K/ J. 


Let us look into a behavior of the lower band in more 
detail. We first obtain extrema via = 0, which 

leads to 


sin k 




= 0 . 

( 10 ) 


If (/) TT with Ff/J ^ 1, fc = 0 and ±7r are the solution 
of Eq. (101, and k = 0 gives the minimum of the band. 
As in the case of the normal cosine band, the dispersion 
near A: = 0 is approximated to be quadratic in k. It is 
straightforwardly shown that the condition for the single¬ 
minimum structure can be expressed as |22l I23j 


(I) 

v^'^/ 1 - sin^ y^ 


( 11 ) 


As K/J decreases, on the other hand, the double-well 
structure starts to show up. The critical point between 
the single- and double-minimum structure is given by 


JkV sin^f) 

V2jj "i_sinyy’ 


( 12 ) 


which means the equality limit of the inequality in 
Eq. (11). Then the dispersion near the bottom becomes 
quartic in k. In the double-well structure case, two sepa¬ 
rated minima ±Q can be obtained from the factor of the 
square bracket in Eq. (fTol) [5^ [23] , 


Q = sin ^ 



In addition, the dispersions around the minimum point 
are quadratic in k as in the case of the single minimum. 
We emphasize that ^ = tt is special because the sym¬ 
metric double-well structure is strongly protected and its 
minima are located at Q = regardless of K/J. 

The similar behavior is also found in the case of 
changes of (j) with a fixed K/J. In that case it is shown 
that for the small enough (j), the band forms single¬ 
minimum structure, and shows up the double-well struc¬ 
ture when (j) goes through a critical value (j)c- A crit¬ 
ical flux (j)c between these two structures is shown to 
be [22l|23| 



Based on the change of the band structure discussed 
above, let us incorporate interaction effects within a 
weakly-interacting regime, J, K ^ U. In this interaction 
regime, the bottoms of the band are also important for 
bosonic systems, which is an essentially different point 
from the fermion systems. Therefore, we first truncate 
all the effects involving the upper band (afc). The Hamil¬ 
tonian ([^ is reduced to 

k ki,k2,k2,k4 

(15) 


where 

^ki,k2,ks,k4 


^fcl+/C2—^3 —^4,27771' 

n'ez 



(16) 


Note that in contrast with a system in continuum space, 
we need to consider scattering processes involving a finite 
momentum transfer equal to the reciprocal lattice vector. 
In our model, the finite momentum transfer to 27rn' with 
an integer n' is allowed, which is nothing but the umklapp 
process, and turns out to play a crucial role in the (j) = tt 
case. 

In what follows, we separately look into many-body 
ground states in each topology of the single-particle band. 


III. BAND WITH A SINGLE MINIMUM 

For weakly-interacting bosons in higher dimensions, 
the Gross-Pitaevskii (GP) approach as one of the mean- 
field theories is known to provide good results [17]. As a 
consequence, the system undergoes a Bose-Einstein con¬ 
densate (BEG), which also implies spontaneous breaking 
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of U{1) symmetry. Thereby, a gapless excitation mode 
known as a Nambu-Goldstone (NG) mode is obtained. 

However, it is also well known that for an interact¬ 
ing one-dimensional bosonic system, there exists nei¬ 
ther BEG nor NG mode in the thermodynamic limit. 
Namely the mean-field analysis underestimates fluctua¬ 
tion effects, and cannot thus capture correct results in 
the one-dimensional cases. However, one-dimensional su¬ 
perfluids in the weakly-interacting regime show very slow 
power-law decay in such a way that the order effect works 
almost comparably to the quantum fluctuation. In addi¬ 
tion, there is also an acoustic phonon mode similar to the 
NG mode, although it does not correspond to the spon¬ 
taneous symmetry breaking. From these facts, it is found 
that the GP approach is not correct in a strict sense, but 
would provide a practically reasonable starting point to 
discuss the ground state and low-energy excitation struc¬ 
ture. In addition the advantage of the GP approach is 
that both kinetic and interaction energies can be simul¬ 
taneously taken into account at the mean-field level 

Based on the above observations, let us consider the 
system with the single band. m As far as the weakly- 
interacting bosons are concerned, the bosons dominantly 
populate the minimum of the lower energy band, and one 
can perform an approximation such that all the energy 
states except for ones in the vicinity of the minimum 
are projected out. Thus the low-energy single-particle 
spectrum is approximated as 

E.{k)^-Eo + ^, (17) 


where Eo = K 


2Jcos and ^ Since 


all the wave numbers are restricted to be \kj\ 1 due 
to the long-wave-length approximation, the effective in¬ 
teraction parameter is approxi¬ 

mated in the following way: The small wave number 
k leads to sin(^fe/2) « and cos(^fc/2) « l/-\/2, 

and by substituting the approximated form of into 
Eq. (16), the interaction parameter is approximated as 
~ U/2. Thus, we reach the following low- 
energy effective Hamiltonian (15): 


H 


dx 


- -f {x)l3{x)p\x)P{x) 


(18) 


where P{x) = ■ We note that this is essen¬ 

tially identical to the Lieb-Liniger model [3^. As far as 
the weak-coupling limit is concerned, one may consider 
the following GP ground state: 


\GS) = ^{PUf\Q), ( 19 ) 

where N is the number of bosons. 

Let us next incorporate long-wavelength fluctuations 
which play crucial roles in low-energy properties. To this 


end, we adopt the hydrodynamic approach also known 
as the bosonization for bosons [26ll28) : 


P{x) 


Vvj(a;) 


E' 

m^’L 


^2im[-Knx— (p{x)\ —i9{x) 


,( 20 ) 


where n = N/L is the mean density. We introduced 
the density and phase fluctuations, ip{x) and 0{x), re¬ 
spectively, and the commutation relation between them 
is given by 


[0{x), = i5{x — x'). 


( 21 ) 


Applying the above bosonization formula (20) to 
Eq. ([l^, we obtain 


~ 2I 


Ao 


( 22 ) 


where Vq = y and Kq = This is the 

Hamiltonian for the celebrated Tomonaga-Luttinger liq¬ 
uid (TLL), which corresponds to a c = 1 conformal field 
theory. It is remarkable that in this TLL Hamiltonian, 
the long-range order (LRO) of the single-particle den¬ 
sity matrix, incorrectly predicted by the GP mean-field 
theory, is directly confirmed to be modified into a cor¬ 
rect quasi-LRO of the algebraic decay ; {P^ {x )pm ~ 
(y)l/2/Co^ 

Let us next look into the rung (^ and chiral current (|^ 
by translating them in effective theory derived above. By 
using the bosonization formula ( |20| the current opera¬ 
tors, Eqs. (§ and are expressed as 


j^{x) ~ 0, 

jc[x) -- 2nJ sin + 0(V^6»), 


(23) 

(24) 


where we point out that the rung current vanishes re¬ 
gardless of bosonization, while the chiral current has 
the nonzero constant term and the terms starting from 
Note that at the level of the long-wave approxi¬ 
mation, a fluctuation of the rung current disappears. 
On the other hand, a fluctuation of the chiral current, 
5jc = jc — ‘2-nJ sin ((()/2) ~ V^0, behaves as 

{5jc{,x)5jc{P)) ~ l/x'^, (25) 


which is given by the Gaussian property in the TLL 
Hamiltonian such that (V^0(a:)V^d(O)) ^ 1/x^. In addi¬ 
tion, Eq. (24) shows that the chiral current increases with 
(j). These properties correspond to the Meissner phase 
introduced in Ref. [S], which has been derived in a con¬ 
dition J ^ K,U different from the present case. 


IV. BAND WITH DOUBLE MINIMA 

We next examine low-energy properties of the system 
with the double-well band structure where we can distin¬ 
guish a commensurate wave number Q, giving the lowest- 
energy single-particle states, from incommensurate one. 
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In the commensurate cases, Q can be represented as 
Q = TTp/q^ where p, q are coprime numbers. 

The effect of commensurability is related to types of 
interactions. As pointed out in Ref. [S], a g-body inter¬ 
action produces the serious effect for Q = np/q. Thus, 
if arbitrary multi-body interactions come into the low- 
energy effective theory, every commensurability should 
be taken care. Note that it does not mean that multi¬ 
body interactions are required at the microscopic level. 
Namely even if only two-body interactions are assumed in 
the microscopic Hamiltonian, multi-body interactions are 
generated as virtual multiple-scattering processes when 
we integrate out irrelevant high-energy degrees of free¬ 
dom such as deriving a low-energy effective Hamiltonian 
and implementing a perturbative renormalization group 
theory. 

However, such virtual multiple-scattering processes 
would be suppressed in the weakly-interacting case, and 
the relevant case would be only for Q = 7r/2 in which the 
two-body interaction yields the strong commensurability 
effect. 

Here we first consider an incommensurate Q case in 
Sec. EXi Next in Sec. EYl we move on to the discus¬ 
sion of the Q = ^ {(j) = tt) case as one of the commen¬ 
surate cases. The other commensurability is also briefly 
discussed in Sec. IIV Cl 


A. Incommensurate Q case 


In contrast with the single minimum case, the mean 
ground-state density with the double well structure de¬ 
pends on the couplings of the Hamiltonian. Following the 
analysis for a BEC on a double-well potential [m , we as¬ 
sume the following ansatz, first introduced in Ref. |16) : 


\GS) = cos7/?^ -b e*®- sin7^lQ)^|0),(26) 


where 7 and 0± are variational parameters. 

By taking the expectation value of H in Eq. 
the above ansatz, one obtains [16] 


(15) with 


Eoh,0±) 

N 


E_{Q) + 


Un 



sin^ 27 


{n+) 

(n_) 


- sin^ + 2 

5 

(27) 

{0l[x)j5+{x)) = ncos^7, 

(28) 

(/3l(a;)/?_(a;)) =nsin^7, 

(29) 


where j3±{x) = ■ We note that these 

mean-field values have no dependence in which im¬ 
plies that the problem is reduced to optimization of the 
single variational parameter, 7. Then, the optimized 7 is 
alternatively determined by whether [16] 


or 


sin^^Q > 1- 


(31) 


In the former case (301, the ground state is minimized 
by 7 = 7r/4 [TS|, where the populations at fc = ±Q are 
the same: (n+) = (u._). Thus, at the mean-field level, 
we expect that there are two independent BECs in the 
ground state. We note that this is different from a BEC 
on a double-well potential, where the ground-state energy 
depends on the relative phase via a hopping term between 
the condensates m- By contrast, such a hopping does 
not exist in our system, and thus the mean-field ground 
state is free to the relative phase. We will come back 
to this point in an analysis in the </) = tt case, where 
a relative phase dependence shows up via the umklapp 
process in a nontrivial manner. 

In the latter cas e (|M| ), the ground state is characterized 
by 7 = 0 or 7r/2 |16|. where the mean density becomes 
((n+), (u-)) = (u, 0) or (0,n). This is the solution such 
that all the bosons occupy either at fc = Q or at A: = —Q. 
Thus, the mean-field theory shows that Z 2 symmetry is 
spontaneously broken in the ground state, and a single 
BEC occurs simultaneously. 

The transition between these mean-field ground states 
occurs at | sin^.jQ = 1, which turns out to be rewritten 
as 


-)■ 

2J 


sin'*((^/2) 


(32) 


What is important is that the above critical KjJ is 
smaller than another critical K/ J, given in Eq. @, be¬ 
tween the single- and double-minimum band topology. 
Namely it means that the solution (32) always exists in 


the regime of A'/2 J and (j) in which the double-minimum 
band structure comes out. Therefore, we see that the 
transition between the mean-field ground states always 
occurs at a certain K/2J given by c/i > 0. 

Let us next look into fluctuation effects based on the 
above mean-field analyses. As in the case of the single¬ 
minimum band, we approximate the Hamiltonian as |I3j 


H ■. 


dx 




17(2-1- sin^ ^q) 


(n+ -b n-Y + 


U{2-isii?^Q) 


(n+-n-f 
(33) 


where effective mass. 

We first incorporate fluctuation effects in the case of 
the mean density (n+) = {'n+) = n/2 which is stable 
when Eq. (30) is obeyed. By using bosonization for¬ 
mula (|20|) in f3±, 


P±{x) 


n 

12 


Vip±{x) 


^ \ ^ 2im[Knx/2 — ip±{x)\ —i9±{x) 


^sin^^Q < 1, 


(30) 


(34) 
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with the density and phase fluctuation fields in the vicin¬ 
ity of the bottoms k = iQ, Lp± and we obtain 


Hes= 


2tt 


i/'=s,a 


dx 




(V^i 

K, \ 


2g 

(27rQ!)2 


dxCOs{ViLpa), (35) 


where we have introduced the symmetric (anti¬ 
symmetric) fields, ifs{a) = ^(‘P-t + (-)‘/?-) and e^^a) = 

^(0+ -I- (—)0_). For convenience’s sake, we have also 
introduced the cutoff parameter a = l/(7rn) H?]- The 
velocities, Vs and Ua, and TLL parameters, and iFo, 
appearing in the quadratic parts of the Hamiltonian are 


specified by 


Va = 
Ks = 
K„ = 


I nil{2 + sin^ ^q) 

AM* ’ 

/ nil{2 — 3sin^ ^q) 
AM* ’ 
n 

M*U{2 + sm'^ ^q)’ 
ri 

M*U{2-3sm^ ^q)’ 


(36) 

(37) 

(38) 

(39) 


and the coupling of the cosine term is given by 

g = Usm'^^Q. (40) 


In the effective Hamiltonian (35), the symmetric and 
anti-symmetric fields are decoupled. Especially, the sym¬ 
metric part is the conventional TLL Hamiltonian. On 
the other hand, the anti-symmetric part seems not to be 
the TLL Hamiltonian due to the presence of the cosine 
term. Thus, the low-energy properties of the system are 
determined by the relevancy of this cosine term in the 
sense of the renormalization group. To this end, we im¬ 
plement the perturbative renormalization group analysis. 
By treating the coupling constant g as a perturbative pa¬ 
rameter, we obtain (See Appendix.), 


^i^=2{l-Ka)g/va, (41) 

where I is the scaling parameter. In general, ATg, Ka >> 1 
in weakly-coupling bosons with a contact interaction. 
This means that g is the irrelevant coupling, and the 
cosine term goes away in the low-energy limit. There¬ 
fore, this phase is found to be characterized by the two- 
independent TLL. 

Let us next look into the currents. The bosonized ex¬ 
pressions of the operators are summarized as 


sin 


jc(x) - nJ 4sin^ ( 

+V2cos^ ( ^ 


j^{x) ~ nK ( sin^ ( %■ ) ~ cos 


-I- Q ) — sin^ 


cos 


Q]^0a 


cos I - - 


V6>a - 4sin cos sin cos (q{2x -f 1) - V29a^ 
sin(2Qa; - V29a). 


(42) 

(43) 


By taking the averages of the these quantities, we obtain 
Ucix)) ~ 4nJsin^ sin (44) 

(j-^(a:))~0, (45) 


where we used (V0a) 
(sin(2Qa; — \/29a) = 0. 


= (cos(Q(2a; + 1) - '/29a)) = 
The form of Eqs. (44) and (45) 


is the same as what Wei and Mueller [T^ have derived 
within the mean-field analysis for the net chiral cur¬ 
rent [29]. As shown in Refs. 0[I6], the chiral current 
monotonically decreases with (j), and goes to zero as ap¬ 
proaching 0 —>■ TT. This phase corresponds to the (in¬ 
commensurate) vortex phase first introduced in Ref. [5] 
in which the reduction of the chiral current has been at- 






















7 


tributed to the penetration of the vortices. Indeed, a sig- current fluctuations, 5jc = jc — (jc) and 6j-^ = 

nature of the vortices is found in the correlations of the They are calculated as 


((5jc(a;)(5jc(0)) 


2 t 2 r 


nV 

~Kn 


sin^ f ^ j ( ^ + <3 


+in^j‘^s\T? ( ^ j cos^ j sin^ j cos(2Qa:)-^^^^_^ 


cos^ f ^ j Q 


1 


2 1 




2 I \ 2 

Sin 1 ^ I — cos 


1 2 


cos(2(5a 


rl/K^ 


(46) 

(47) 


where we used the facts that (V6*(a;)V0(O)) ~ 
^gAi6»(a;)g-Aie(o)^ ^ ^ Constant A, and the 

X 

cross terms such as {\/9{x) cos{Q — •\/20a(O))) vanish. 
They have now oscillation components decaying with a 
power law. Recalling Ka ^ 1, these power-law decays 
are extremely slow. 

We next consider fluctuation effects in the case of the 
biased mean density ((n+), («-)) = (n, 0) or (0,n). 
For the sake of simplicity, let us take the case of 
((n_|_), (n_)) = (n,0) |30]- Namely all the bosons only 
populate around k = Q. The degrees of freedom around 
k = —Q are completely suppressed as far as the inter¬ 
action is weak enough, and we may consider only the 
degrees of freedom around k = Q. Then the Hamilto¬ 
nian simplified by the long-wave-length approximation is 
given as 


H = 



Y 72 

/3^ {x)———j5(.x) 
^ 2M* ^ 


f/(2-sin^eQ) ^2 


48) 


Furthermore, by the bosonization formula (20), we obtain 

\K{vef 


i7eff=^/ dx 


^2, 


K 


(49) 


where the velocity and TLL parameter are, respectively. 


V = 


K = 


I nil {2 — sin^ ^g) 

TM* ’ 

/ 2n 

/17*17(2-sin^eJI' 


(50) 

(51) 


Unlike the equal density case, the effective theory is de¬ 
scribed by a single TLL. 

Let us next look at the currents. They are bosonized 
as given by 


jc{x) ^ AnJ 



y{x) - 0, 


— 2nJ 





Ve, (52) 


(53) 


where the rung current is shown to be zero at the level of long-wave approximation as in the case of the Meissner 
phase. The averages of them are calculated as 


{jc{x)) - 4nJsin^ sin + <3^ > 

(j-^(a:)) - 0, 


(54) 

(55) 


which are identical to the expressions for the net currents obtained in Ref. [T^ and look the same as those of the 
(n+) = (n_) case, i.e., Eqs. ([44|) and (45). However, this does not mean that all the low-energy properties coincide 


with the equal-density case. Indeed, we find that the difference occurs in the current fluctuations as 


2?7^ 

{d3c(x)53M) - 
{5y{x)5y{Q)) ~ 0. 


sin^ 





(56) 

(57) 


Thus, in contrast with the vortex phase in the (n+) 


(n_) case, the current fluctuations in this phase do not 

































\m\ln 



KI2J 


FIG. 3. The absolute value of density difference at ^ 
in the biased ladder phase. At the boundary between biased 
ladder (m 7 ^ 0) and Meissner phases (m = 0), the density 
difference disappears. On the other hand, at the boundary 
between incommensurate vortex and biased ladder phases, an 
infinite number of degeneracy in density difference emerges 
due to the emergent symmetry as shown in Sec. |IVA1| 


have an oscillating component. 

A peculiarity of this phase is seen in the density in 
each leg, ni and 712 . To see this, we define a magne¬ 
tization, i.e., population imbalance between the legs, as 
m = (ni) —( 712 ). According to the mean-field theory (26), 
the magnetization can be calculated as m 


= -71 cos ^Q. 


(58) 


We note that it can be proved that the magnetization 
value m is robust even when the quantum fluctuation up 
to the bosonization level is incorporated. This is due to 
the fact that the fluctuation of the magnetization Sm is 
given as Sm ~ (Vtp) = 0. Thus, the mean-field result 
is applicable. As shown in Fig. m takes a nonzero 
value as far as the stability condition of the phase is met, 
which means the spontaneous imbalance of the popula¬ 
tions between the legs occurs. This phase corresponds to 
the biased ladder phase introduced in Ref. m, which has 
been first demonstrated within the GP mean-field theory. 
What is addressed here is that even in the presence of the 
quantum fluctuation the Z 2 symmetry breaking between 
the populations in the doubly-fold lowest energy in the 
double-well band is maintained, and the biased ladder 
phase is thus stable at the full quantum level. This would 
be reasonable once one recalls the fact that even though 
a continuous 17(1) symmetry cannot be broken in quan¬ 
tum one-dimensional systems, a spontaneous breaking of 
a discrete symmetry is possible. 


1. Analogy with ferromagnetic XXZ model 


Let us now examine the nature in the phases and tran¬ 
sitions between them, from the viewpoint of symmetry. 
To this end, we focus on the low-energy Hamiltonian (331, 


which is invariant under the continuous transformations, 
/3± —>■ and discrete transformation, (3± —)• 

This implies that symmetry of the low-energy Hamilto¬ 
nian (33) is 17(1)+ X 17(1)_ X Z 2 , where the subscript 
± of 17(1) represents the corresponding symmetry in /3+. 
We also note that the above symmetry can also be rep¬ 
resented as 17(1)V X 17(1)A X Z2 where 17(1)+ and 17(1) a 
represent the vector 17(1) symmetry, f3± —>■ e'^^l3± 


■ and 


the axial 17(1) symmetry, P± —>• respectively [HT] . 

In the vortex phase, the low-energy effective properties 
are captured by the two independent TLLs, in which the 
two independent 17(1) symmetries are hold; namely no 
symmetry breaking occurs in this phase. In the biased 
ladder phase, on the other hand, the low-energy effec¬ 
tive properties are described by the single TLL reflecting 
the acoustic phonon excitation around one of the two 
minimum-energy states in the double-well band. The Z 2 
symmetry turns out to be broken, since the minima of 
the band are degenerate and one of the minima is spon¬ 
taneously chosen. 

It is important to clarify the transition point between 
the vortex and biased ladder phases, where we need a 
special attention. First, we address that at this transi¬ 
tion point, symmetry of the low-energy Hamiltonian is 
enlarged. To see this, it is convenient to introduce a two- 
component spinor. 


/3 = 


/ 3 + 

/3- 


(59) 


In this spinor representation, we can define 17(1)+ x 
SU{2) transformations as /3 —>■ 

where ai (i = x, y, z) is the Pauli matrix, and ai {i = 
1, • • • ,4) represents an angle of 17(1)+ x SU{2). Then, 
the low-energy Hamiltonian (33) is shown to be invari¬ 
ant under the above transformations at the transition 
point where the term proportional to (77+ — 77 _)^ disap¬ 
pears, i.e., the Hamiltonian has 17(l)+xS'17(2) symmetry. 
This implies that we can define the following conserved 
charges: 


S+ = J dx(3ll3., 
S-= j dxl3lp+, 

Sz = J dx[n+ — 7i_], 
Ny = J dx[n++n-]. 


(60) 

(61) 

(62) 

(63) 


Here, Eqs. (60)-(62) constitute SU{2) charges and 


Eq. (63) originates from 17(1)+ symmetry. In addition. 


the Hamiltonian (33) at the transition point is identical 


to the two-component bosonic Yang-Gaudin model where 
the Bethe ansatz solution is available [5^1 - 151] . So far 
the followings on two-component bosonic Yang-Gaudin 
model are known: As usual, there is a TLL associated 
with 17(1)+ symmetry. In addition, the SU{2) symme¬ 


try corresponding to Eqs. (60)-(62) is spontaneously bro¬ 
ken in the ground state, which is essentially corresponds 
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to the physics of the Heisenberg ferromagnet 
Due to the spontaneous breaking of the SU{2) symme¬ 
try, we expect that there is a NG mode whose disper¬ 
sion is quadratic in k and there exist an infinite number 
of the degenerate ground states and one of them is se¬ 
lected spontaneously. This is significant difference from 
the biased ladder phase where there only exists double 
degeneracy. 

The scenario discussed above reminds us of the similar¬ 
ity to the ferromagnetic XXZ model. When we bosonize 
the XXZ model, the low-energy effective theory is de¬ 
scribed as the sine-Gordon model. When the XY in¬ 
plane anisotropy is strong, the theory is renormalized to 
the TLL which corresponds to the so-called XY phase, 
and upon approaching the isotropic point the velocity 
and Luttinger parameter of the XY phase, respectively, 
go to zero and infinity, which implies the quadratic dis¬ 
persion of an excitation |27j . Furthermore, beyond the 
isotropic point, i.e. Ising anisotropy, the XXZ model 
undergoes the so-called ferromagnetic Ising phase, where 
the excitations are massive, and the Z 2 spin-inverse sym¬ 
metry is spontaneously broken. 

In our case, the anti-symmetric sector in the vor¬ 
tex phase (351 exhibits the same effective theory as the 
XXZ model, and the velocity (37) and Luttinger param¬ 
eter (39) show the same behavior as those in the isotropic 
limit of the ferromagnetic XXZ model. Going beyond 
the Heisenberg point, we encounter the biased ladder 
phase described by a single TLL, which corresponds to 
the ferromagnetic Ising phase in the ferromagnetic XXZ 
model. It is interpreted that for the biased ladder phase 
the anti-symmetric sector goes away into the high enegy 
regime, which would correspond to the massive excitation 
in the ferromagnetic Ising phase. Namely such a massive 
excitation should describe the change of the populations 
on the band bottoms and is regarded as the high energy 
one in our approach, which is excluded in Eq. (49). In 
Sec. IV A 2[ it is shown that such a massive excitation 
may be incorporated with the Bogoliubov theory. 


2. Bogoliubov spectrum in biased ladder phase 

To see some insight into the biased ladder phase, let us 
here review the Bogoliubov theory given by m- Since 
the Bogoliubov theory is based on the expansion from the 
GP solution, it must underestimate fluctuation effects. 
However, this does not mean that all the results by means 
of the Bogoliubov theory are incorrect as pointed out in 
Sec. |HI| At least, an excitation spectral feature in the 
low-energy limit for weakly interacting one-dimensional 
bosons is expected to reproduce the correct behavior. In¬ 
deed, it is known that linear excitation feature occur¬ 
ring in the Lieb-Liniger model [35] and quadratic ex¬ 
citation feature occurring in the two-component Yang- 
Gaudin model [33] can be captured by the Bogoliubov 
theory. 

As usual, by applying Pk = \^Sk,Q + Pk, where Pk is 


the fluctuation field and Nq is the number of the particles 
in the condensate, to Eq. (15), we obtain the Bogoliubov 
Hamiltonian to be correct up to second order in the fluc¬ 
tuation field 




fc >0 


— (^Q+fcJ ^Q-k) 


Cik) rj{k) 
r]{k) ((-k) 


I^Q+k 

'g-fc 


(64) 


where the matrix elements are defined as 


C(fc) = E_{Q + k)- E_ (Q) -h 2Un 


• 2 -2 ^Q+k 

sm — sm —^— 
2 2 


2 2 ^Q+k 

- cos — cos — 
2 2 


— Un 


sm 


ri{k) = Un 


■ 2 ■ ^Q+k . 

sm — sm —^— sm 
2 2 


2 

ig-k 


cos^ , (65) 


I 2 ^Q+k ^Q—k 

+ COS — cos -cos - 

2 2 2 

Here, we did an approximation N 
ize the above Hamiltonian (64), a simple scheme is to 
consider the following eigenvalue problem |351136j : 


( 66 ) 

Nq. To diagonal- 


1 0 \ /C(fc) T]{k) 

0 - 1 ] U(fc) c(-fc) 


(67) 


Then, the Bogoliubov spectra ±e are obtained as eigen¬ 
values of the above matrix equation. The Bogoliubov 
spectrum is shown to be m 


( 68 ) 

A behavior of the Bogoliubov spectrum is shown in Fig.[^ 
When looking at fc —)■ 0, we have the linear spectrum, 
which can be regarded as the TLL. In addition, there is a 
local minimum around k = —Q. This roton-like behavior 
would be interpreted as the massive excitation originat¬ 
ing from the Z 2 symmetry breaking. Indeed the similar 
situation recently realized in Ref. m also occurs in a 
BEG in a shaken optical lattice, in which Z 2 symmetry 
is broken. The energy gap around k « —Q, e{k « —Q), 
is found to go closed as approaching the point where the 
symmetry is enlarged to SU(2). This scenario may re¬ 
mind one of the ferromagnetic transition in the XXZ 
model when approaching the Heisenberg point from the 
Ising anisotropic side, the Ising gap will collapse and the 
dispersion turns to quadratic in k. 


B. Commensurate Q case {Q — V2) 

Let us next consider a commensurability effect, i.e. 
Q = 7r/2 ((/) = tt) case where we require considerable 
attention. In this case, the double-well band struc¬ 
ture is always maintained and its minima are located 
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e(k) 



termined by 7 = 7 r/ 4 , which is independent oiKjJ. This 
ground state leads to the balanced density (n+) = (n_). 
Thus, the biased ladder phase is washed out in the pres¬ 
ence of the umklapp process. 

Let us next consider quantum fluctuations from the 
mean-field solution. In this case, we need to retain the 


following process in the effective Hamiltonian (331: 


FIG. 4. Schematic behavior of the Bogoliubov spectrum pre¬ 
dicted by Eq. ( | 68 [ ). The excitation is linear in small fc, which 
should be interpreted as the TLL spectrum. On the other 
hand, a local minimum exists in the vicinity of fc = —Q orig¬ 
inating from the Z 2 symmetry breaking. 


a,t k = ± 7 r /2 regardless of the ratio K/J. A peculiarity 
is an emergence of the umklapp scattering process be¬ 
tween two energy minima in the band, which has been 
overlooked so far. To consider this, we turn back to the 


iJumkiapp = u sin^^Q J dx[/3l_l3l_/3-f3- + h.c.]. (71) 

Note that due to this term symmetry of the Hamil¬ 
tonian is lowered from C/(l) x 1 /( 1 ) x Z 2 to U{l)v x ^ 2 - 
Thus, the axial 1/(1) symmetry {(3± —>■ e^^^(3±) disap¬ 
pears from the low-energy Hamiltonian, and the contin¬ 
uous symmetry remaining turns out to be only the vector 
17(1) symmetry (/3± —>■ e^^(3±). On the other hand, the 
Z 2 symmetry {j3± —> jdz^) remains in the presence of the 
umklapp term (In])- 

Let us next perform bosonization as follows: 

^ = E S / dx\K^{we^f ' 


K„ 


^—s,a 

91 I 


(27rQ!)^ 


dx cos{V8(pa) — 


92 


(2t 


dxCOS^VSOa), 

(72) 


mean-field ansatz (26). The form of this ansatz includes where gi — sin ^q/2, g 2 — sin and 


the couplings of the far-separated states k = ±Q, and 
thus automatically allows us to take into account the 
umklapp process involving the large momentum trans¬ 
fer. Based on the mean-field ansatz, the ground-state 
energy is straightforwardly calculated as 


EQ{-f,e±) _ Un 

3 , cos(26»+-26i_) 
2 2 


M* 


+ 


- sin^ Cq + 2 
sin^ ~ 11 sin^ 27 .(69) 


K, = 


Ka = 


n 


4M*U 


Note that it differs from Eq. (27) due to the presence of 
the umklap scattering. As can be clear from the above 
expression, the energy must be minimized when the rel¬ 
ative phase satisfies 0+ — = ± f. We notice that this 

is different from the case of a BEG in a double-well po¬ 
tential, where the relative phase is zero in the ground 
state m- This difference originates from the fact that 
the relative phase dependence is caused by a hopping 
(kinetic) term such as — Jcos(0+ — 0_) in a BEG on a 
double-well potential while that in our model originates 
from the interaction term in our model. 

By substituting this relative phase 0+ — 0_ = 7 r /2 or 
— 7 r /2 into Eq. (69), the ground-state energy is going to 
be 


4M*C/(l-sin^^Q) 


(73) 

(74) 

(75) 

(76) 


An essential difference from the incommensurate cj) case 
is the presence of cos VSda which comes from Tlumkiapp ■ 
Eurthermore we move on to the renormalization group 
analysis to see the low-energy properties of the system. 
The parameters in Eq. (72) are found to obey the follow¬ 


ing renormalization group equations as (See Appendix) 
d{gi/va 


dl 

d{g2/va) 

dl 


= 2(1 - Ka)gi/Va, (77) 

= 2(l-^)g2/va. (78) 


^^ 0 ( 7 ) 


N =^-«5) + Xl 


(sin^ - 1 ) sin^ 27 


— Sin 


■Cq + 2 


(70) 


Since Ka 3 > 1 in the weakly interacting case assumed 
here, it is found from Eqs. (77) and (78) that gi/va 
and g 2 lva are rapidly renormalized, respectively, to be¬ 
ing zero and divergent as I increases. Namely, cos \/80a 
in Eq. (72) is highly relevantly retained in the effective 


Since sin^ < 1, the ground state can be uniquely de- 


Hamiltonian in the low-energy limit while cos V8(t>a ir¬ 
relevantly goes away as in the case of Q 7 ^ tt/2: The 
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(a) 

Incommensurate 

vortex 


Biased ladder 


Meissner 


KIJ 

(b) 

Commensurate vortex 


KIJ 

FIG. 5. Schematic phase diagram (a) for Q ^ 7r/2 and (b) 
for Q = 7r/2. An emergent SU{2) symmetry shows up at the 
boundary between incommensurate vortex and biased ladder 
phases, which is marked with the red circle in (a). 


renormalized effective theory is 




fi—s,a 


'J± 

27T 


dx 




■( 2 ' 


[ dxcos{V89a). (79) 


Thus, the anti-symmetric sector becomes gapful due to 
the fixed relative phase (cos(-\/80a)) = 1, which is consis¬ 
tent with the mean-field scenario, deduced by Eq. (691, 
with the umklapp scattering, 

Let us also look at the currents in this phase. By means 
of the bosonization formula, the averaged currents are 
calculated as m 


{jc{x)) ^ 2nJ sin^Q (-1)^, (80) 

{j^ix))nK cos (81) 


Thus, the current pattern shows up due to the umklapp 
effects. In a similar manner, we can evaluate the current 
correlations, which turn out to be zero (or exponentially 
decay in x). This phase is called (commensurate) vortex 
or chiral superfluid phase. 

The possible phase diagram in the weak coupling 
regime J, K ^ U is summarized in Fig. 


1. Transition between commensurate and incommensulate 
fluxes 


Now we discuss properties in the vicinity of the com¬ 
mensurate flux (f> = TT. To this end, we first estimate 
the gap coming from the coupling cos-\/80q in Eq. (79) 


based on the renormalization group equation |27j . As 
mentioned above, this coupling g 2 is highly relevant in 
Ka 3> 1 and is therefore expected to flow to the strong 
coupling regime (52 —t 00 ). Thus, we introduce a typical 
length scale I* estimated from Eq. (78) as 


( 82 ) 


and stop the renormalization group flow at the length 
scale I*. On the other hand, if we are in the strong 
coupling, we may approximate the cosine term as 
52 cos(-\/80a) ~ 52(1 — ^(^a{x)]- Therefore, applying the 
expansion to the renormalized effective theory (79), the 
Hamiltonian is easily diagonalized, and the consequently 
estimated gap A(Z*) for the renormalized coupling con¬ 
stant 52 (^*) at a termination of the renormalization group 

flow is found to be A{1*) ^ ^ \J Taking 

into account that the gap is renormalized as A(Z) = e*A, 
we obtain 



Based on this estimation, we next consider the situa¬ 
tion in which Q is slightly deviated from a commensurate 
point Q = 7r/2, and write Q as Q = f + f j where S is 
assumed to be small. Then, the umklapp term 0 is 
bosonized as 

77umklapp ~ -TTT^ / dx COs{VS0a + Sx), (84) 

(27ra)^ J 

which is namely the oscillating term as a function of 
X. We may safely make the replacement of <5 — ?► 0 if 
S is irrelevant. However, if 6 is relevant, the umklapp 
term experiences a strong oscillation, and cancels out in 
the renormalized Hamiltonian. In general, the transition 
from the irrelevant and relevant S or vice versa is known 
to occur around A ^ <5, and called a commensurate- 
incommensurate transition |27] . In our cases, the transi¬ 
tion between the chiral superfluid and vortex phases cor¬ 
responds to such a commensurate-incommensurate tran¬ 
sition [27], because once <5 becomes relevant, the effec¬ 
tive theory reduces to that of the two independent TLLs, 
which is identical to that of the vortex phase. 

On the other hand, the transition between the com¬ 
mensurate vortex and biased ladder phases may not 
be captured by this scenario of the incommensurate- 
commensurate transition because the effective theory in 
the biased ladder phase is not the two independent TLLs. 

It is also interesting to consider the transition between 
the commensurate vortex and Meissner phases. However, 
it is difficult to describe such a transition by means of our 
approach. Thus, it would be worthwhile examining the 
nature of these transitions in a numerical simulation. 


C. Other commensurability effect 

So far, we have discussed only the (p = tt case as a 
commensurate case. The point on the p = n case is that 
the biased ladder phase is suppressed by the presence of 
the umklapp scattering. Here we consider the other com¬ 
mensurate case in order to see roles of general commensu¬ 
rability, and we will see that the (p = tt commensurability 
is special. 
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As heretofore, let us first consider the mean-field ap¬ 
proximation for the Hamiltonian (15). Then, it turns 
out that as far as such a Hamiltonian is concerned, we 
always have the energy expression (27) regardless of com- 
mensurability of Q. Thus, the mean-field phase diagram 
is identical to the incommensurate Q case. 

Let us next consider the bosonization to see quantum 
fluctuation effects. For the incommensurate vortex phase 
at the mean-field level, it is expected that the higher or¬ 
der perturbation theory generates cosine terms in Oa at 
a commensurate Q as shown in Ref. [9]. Since many 
of such cosine terms are relevant for Ka S> 1, we find 
that the incommensurate vortex phase is replaced by the 
commensurate vortex phase by quantum fluctuation ef¬ 
fects if the coupling of the cosine term is larger than the 
temperature fFf\ . 


On the other hand, for the biased ladder phase at the 
mean-held level, we hnd that a cosine term to hx 6a does 
not show up by the bosonization since the mean density 
in one of the wells is equal to zero. Namely, the biased 
ladder phase at such a commensurate Q is robust. 


V. SUMMARY AND PERSPECTIVE 


We have examined the two-leg Bose-Hubbard ladder 
model subject to a magnetic held hux. We have partic¬ 
ularly revealed the structure of the phase diagram in a 
weak-coupling regime by using a couple of the effective 
theory methods. What we stress is that we have also 
found the so-called biased-ladder phase, hrst predicted 
by the GP mean-held approach, to be robust against 
quantum huctuations. It has also been shown that the 
transition between the biased-ladder and vortex phases 
has a similarity as that of the ferromagnetic XXZ model 
where the emergent SU (2) symmetry comes out at the 
transition between the biased ladder and vortex phases. 
In the case of the ladder system subject to the magnetic 
hux, commensurability works to phase degrees of free¬ 
dom, which produces a kind of umklapp processes. By 
incorporating such an umklapp process at the mean-held 
level, we have shown that the biased-ladder state tends 
to be destabilized by the umklapp process, and turns out 
to be forbidden for the case of </> = tt. 


A. Transition between Meissner and biased ladder 
phases 

As seen in Sec.|H) the dispersion becomes quartic at the 
critical point between the single- and double-minimum 
band structures. In the absence of an interaction, one 
can naively expect that all the bosons condense at the 
lowest energy, and just forms a BEG which is the same 
as the case of the quadratic dispersion. 

A question is what happens in the presence of an in¬ 
teraction. Here we briehy discuss a possible scenario. 


We hrst point out that the similar situation can be also 
considered for one-dimensional two-component bosons 
with spin-orbit couplings, where the bare single-particle 
dispersion becomes quartic at a certain value of the spin- 
orbit coupling and biased chemical potential between the 
two species. By employing the hydrodynamic approach 
and Gaussian approximation, it is shown that the low- 
energy effective theory undergoes non-TLL |35]. Then 
the excitation is still gapless, but is no longer identical to 
an acoustic phonon: the quadratic-dispersion mode. It 
would rather be that of the Heisenberg ferromagnet and 
two-component bosonic Yang-Gaudin model where the 
spontaneous symmetry breaking and NG mode show up. 
Interestingly, however, the off-diagonal density matrix is 
shown to decay exponentially as 

{bl^pbo^p,) ~ (85) 

where is the correlation length given by = 
\J2pq/( ragX^) with a mean-density po, atomic mass to, 
density-density interaction g and spin-orbit coupling A. 
Thus, it means that even one-dimensional superfluidity 
is destroyed. 

From the above example, we can expect the same 
physics in our model. Namely the system might form 
such a non-TLL when the system transits from the Meiss¬ 
ner to the biased ladder state. However, then the Meiss¬ 
ner current would be predicted to be still protected be¬ 
cause its presence is guaranteed by the two band struc¬ 
ture from the ladder geometry and flux (§. m 

In addition to the physical properties, this problem on 
interacting bosons for quartic bare dispersion would have 
another interesting aspect. In general, it is expected that 
an existence of gapless modes supports LRO or quasi- 
LRO while our model is an exceptional case on this state¬ 
ment. Thus, the profound understanding on the role of 
gapless modes in one dimension remains an open ques¬ 
tion. 


B. Stronger U effect 

In the paper, we have analysed the system under the 
condition K, J ^ U. A natural question to come up 
then is what happens when the system with a stronger 
interaction U is concerned. In many of one-dimensional 
systems, both the weak- and strong-coupling analyses are 
continuously connected to each other consequently EH, 
but we would think that the strong coupling regime in¬ 
cludes different physics in our model. This is because 
while in a weak coupling, the low-energy properties can 
be well captured by assuming the quasi-condensates at 
well-separated lowest energy single particle states in the 
double-well band structure, such a picture is no longer 
applicable in a strong coupling due to a hard-core feature 
analogous to Fermi statistics. Namely the double-well 
band feature can be no longer important in low-energy 
physics if we naively assume fermion-like occupation of 
the particles in the band picture as a strong coupling 





13 


limit, and the different properties from those for the weak 
coupling should be then found. Indeed, the recent nu¬ 
merical analysis |39j in a strong coupling shows that the 
biased-ladder phase is not found while the presence of 
the Meissner and vortex phases is confirmed. 

Let us make a further consideration on the physics in 
the intermediate interaction. Then it is convenient to 
take the generalized mean-field ansatz m, 

IGS") = cos 74 -h e*®- siny^) |0),(86) 

where now the wave-vector k pointing at the bottoms 
of the band is also treated as a variational parameter. 
This generalization means that a modification of the 
band structure by an interaction is taken into account. 
As shown in Ref. HU, the variational approach shows 
that the optimized value of k decreases with U, and ap¬ 
proaches zero at a certain Uc- It is not clear whether 
this mean-field ansatz correctly captures the physics in 
the regime U ^ Uc, but we can naively guess, at least. 


from this discussion, that the interaction works so as to 
collapse the double-well band structure. 

Combining the mean-field and numerical result, the 
following scenario can be deduced: The biased-ladder 
state for the weak-interaction goes unstable as U in¬ 
creases; it eventually transits at U = Uc to Meissner 
state, and such a Meissner state continues to that of the 
strong coupling regime which is found in Ref. |39j . To 
test this scenario, or to precisely estimate the critical Uc, 
an unbiased numerical simulation would be necessary. 

Note added: Recently, we noticed a paper m, which 
found the biased ladder phase in a regime J ~ K U 
by means of the density matrix renormalization group. 
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Appendix A: Derivation of renormalization group equation 


In this Appendix, we derive renormalization group equations for the following Hamiltonian: 


H = Ho 

V 


(27rQ!)^ 


dx cos{V8(p) 


(27ra)^ 


dxcos{V89), 


Ho = 


27r 


dx 




(Al) 

(A2) 


where Ai and A 2 are the couplings of the cosine terms, and the cutoff parameter a turns out to play an important 
role in obtaining renormalization group equations. To this end, we adopt a scheme based on correlation functions 
m. which is known to be useful in many one dimensional systems m- 
To be specific, we consider the following correlation function: 

R(ri - r2) = (A3) 

where ri = {xi, y = VTi) {i = I, 2), and the average is taken with the following partition function: 

Z = j (A4) 

If the couplings Ai and A 2 are absent, the Hamiltonian becomes the TLL one (Gaussian), and therefore, we can easily 
evaluate the above correlation function for ri — r 2 ^ a as 




„-ifFi(ri-r2) 


(A5) 


x'^ + {v\T\+a)'^ 


, and ■ ■)ho means that the average is taken with Hq in Eq. (A4). While in the 


where Fi (r) = | In 

presence of the couplings, one cannot evaluate the partition function exactly, at least, one can perform a perturbative 
calculation by assuming that the couplings are small. Then, it is easily to show that the first-order terms in the 
couplings become zero. Thus, up to the second-order terms, the correlation function is given by 

A? 


i?(ri - rz) = 


1-h 


2(27rQ;)^u^ 




e=±l ' 


I ^2 f ,2^/ .2 // -F-iFi(r'-r")/ -2[F2(ri-r')-F2(ri-r")+F2(r2 

2(27rQf)^n^ 


-r")-F2(r2-r')] 


-1) 


(A6) 
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where F 2 {r) = —iArg {vt + a + ix), and cPr = vdxdr. To go further, let us use the fact that dominant contributions 
in the above integrals comes from regions where r' and r" are not too distant. Thus, by introducing the center of 
mass and relative coordinates, R = ’’ and r = r' — r", one can expand the exponential terms in the parentheses 
and obtain 


i?(ri-r2) = 1 

\ 2 [ - R) - Fi(r2 - R)){^\ + V|.)(Fi(ri - R) - Fi(r2 - R)) 

l[Z7Tapv^ J 

- ../L 2 /- R) - F^ir^ - i?))(V| + V|.)(F2(ri - R) - F^ir^ - R)) 

2(2Tra)^v^ J 

By using the so-called Cauchy relations, 

VxFi = lVyF2, 

VyFi = —iVxF2, 


(A7) 

(A8) 

(A9) 


and (V^ -I- Vy) ln(i?) = 2'it5{R), the correlation function is given by 


R[ri - r2) = 


1-f 


>lKFi{ri r2) f 2 2-KF^(r) _ AiAl(ri r2) f 2 2-K-^F^(r) 


This expression allows one to introduce an effective exponent ATefi, 




KeS = K- 


yfK^ r dr f dr /'r\^-K 


Q 


+ ^ - - 


G)' 


. (AlO) 


(All) 


where we introduced j/i = ^ (i = 1, 2). Considering that the effective exponent should not affect by varying the 
cutoff as a' = a -I- da, we obtain 


K{a') = K{a) 
yl[a')=yl{a) 
yl[a')=yl{a) 


y"l{a)K"^{a) da 


2 a 

,, 4 - 4 K{a) 



yl{a) da 
2 a ’ 


(A12) 

(A13) 

(A14) 


By introducing the scaling parameter I, which satisfies a = ooe* with the original cutoff ag, the following renormal¬ 
ization group equations are obtained: 


dK 

~dr 

dyi 

dl 


yl 

2 2 ’ 
2{l-K)y4, 


dy2 

dl 


= 2{1-1/K)y2. 


(A15) 

(A16) 

(A17) 


The low-energy effective theory (35) in the incommensurate vortex phase corresponds to Eq. (41) with the constraint 
A 2 = 0. Thus by fixing 2/2 to be zero in Eqs. (A15)-(A17) and replacing yi = g/va, the renormalization group 
equation (41) is derived. On the other hand, the effective theory (72) in the commensurate vortex phase is identical 
to Eq. (Al), and thus the renormalized group equations, Eqs. (77) and (78), are exactly obtained by the replacement 
yi = giT ^ and y 2 = g^/va- 
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